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ABSTRACT 

Deep images of 10 early-type galaxies in low-density environments have been obtained 
with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The 
global properties of the globular cluster (GC) systems of the galaxies have been derived 
in order to investigate the role of the environment in galaxy formation and evolution. 
Using the ACS Virgo Cluster Survey (ACSVCS) as a high-density counterpart, the 
similarities and differences between the GC properties in high- and low-density en- 
vironments are presented. We find a strong correlation of the GC mean colours and 
the degree of colour bimodality with the host galaxy luminosity in low-density en- 
vironments, in good agreement with high-density environments. In contrast, the GC 
mean colours at a given host luminosity are somewhat bluer {A{g — z) ^ 0.05) than 
those for cluster galaxies, indicating more metal-poor {A[Fe/H] ~ 0.10 — 0.15) and/or 
younger {Aage > 2 Gyr) GC systems than those in dense environments. Furthermore, 
with decreasing host luminosity, the colour bimodality disappears faster, when com- 
pared to galaxies in cluster environments. Our results suggest that: (1) in both high- 
and low-density environments, the mass of the host galaxy has the dominant effect 
on GC system properties, (2) the local environment has only a secondary effect on 
the history of GC system formation, (3) GC formation must be governed by common 
physical processes across a range of environments. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
formation 



1 INTRODUCTION 

The study of extragalactic globular cluster (GC) systems 
is now well-established as a powerful pr obe of the forma- 
tion histories of ellipt ical galaxies (e.g. I West et al.l |2004| : 
iBrodie fc Stradedl2006l l. Because each cluster is essentially 
of a single age, a single metallicity, and a simple stellar 
population (SSP), the colours and line strengths of individ- 
ual clusters are more easily interpreted using evolutionary 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS5-26555. 
t E-mail: sjyoon@galaxy. yonsei. ac.kr 



synthesis stellar population models than is the integrated 
starlight from the spheroidal component. The distribution 
of the aforementioned properties over the entire cluster pop- 
ulation also allows for a unique examination of the relative 
importance of cluster formation episodes as a function of 
age and metallicity. Furthermore, the well-documented cor- 
relations between the properties of GC systems and those of 
their host galaxies (e.g. [Ashman & Zcpf 1998; L arsen et al.l 
'2OOII; iBrodie fc Stradej [2009 : IPeng et all i2006a , hereafter 
PJC06) indicate that a close connection must exist between 
the formation of the diffuse stellar component of a galaxy 
and that of its GC system. 

Most of the extant studies on globular systems in early- 
type galaxies have concentrated on high-density environ- 
ments where these morphological types are most preva- 
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Table 1. The galaxy properties 



Name 


Mb 


RA (J2000) 


DEC (J2000) 


TvDe 


m-M 




E{B -V) 


Density 






(h:m;s) 


(°:':") 










(Mpc^) 


NGC 3818 


-20.33 


11:41:57.36 


-06:09:20.4 


E5 


32.80+0.61 


12.47 


0.036 


0.20 


NGC 7173 


-19.96 


22:02:03.19 


-31:58:25.3 


E+ p 


32.48+0.19 


12.52 


0.026 


0.35 


NGC 1439 


-19.95 


03:44:49.95 


-21:55:14.0 


El 


32.13+0.15 


12.18 


0.030 


0.45 


NGC 1426 


-19.65 


03:42:49.11 


-22:06:30.1 


E4 


31.91+0.18 


12.26 


0.017 


0.66 


NGC 3377 


-19.18 


10:47:42.40 


+13:59:08.3 


E5 


30.25+0.09 


11.07 


0.034 


0.49 


NGC 4033 


-19.11 


12:00:34.74 


-17:50:33.4 


E6 


31.64+0.23 


12.53 


0.047 


0.38 


NGC 1172 


-19.10 


03:01:36.05 


-14:50:11.7 


E+ 


31.66+0.20 


12.56 


0.064 


0.28 


NGC 3156 


-18.84 


10:12:41.25 


+03:07:45.7 


SO 


31.75+0.14 


12.91 


0.034 


0.20 


NGC 3073 


-18.78 


10:00:52.08 


+55:37:07.8 


SABO 


32.64+0.93 


13.86 


0.010 


0.28 


IC 2035 


-18.57 


04:09:01.87 


-45:31:03.1 


SO 


31.18+0.15 


12.61 


0.013 


0.16 



Col. 1: Galaxy name; Col. 2: Absolute B magnitude; Col. 3: RA; Col. 4: Dec; Col. 5: Morphological 
type; Col. 6: Distance modulus; Col. 7: Total apparent B magnitude; Col. 8: Galactic reddening; 
Col. 9: Local density. The morphological t ype and apparent B magnitude are taken from NED. 
The distance moduli are SBF values from iTonrv et all |20o3) except for IC 2035 where the dis- 
tance modulus was estimated using the redshift from NED. The Galactic extinction is ta ken from 
ISchlegel. Finkbeiner. fc David l|l99^ ') and the local density was obtained from ITuUvI l|l98d ). 



lent. Such stud i es in clude the ACS Virgo Cluster Sur- 
yev jCote et al.l |2004) . the ACS Fornax Cluster Survey 
Jordan et al. l2007al ). and the ACS Coma Cluster Survey 



Carter et al.l |2008|). Semi-analytic models of galax y for- 
mation (e.g. iBaugh et al] 11998 : iBeaslev et al] I2OO2I ) pre- 
dict that the largest spread in galaxy properties such as 
metallicity and age should occur for early-type galaxies (L 
< Lt) in low-density environments outside of rich clus- 
ters. Low-luminosity ellipticals also provide a critical test 
for current theories on the origin of bimodal colour dis- 
tributions in GC systems. In some scenarios, clusters in 
the blue (low-metallicity) peak are primarily associated 
with the progenitor galaxy; r ed cluster populat ions form 
late r via multi-phase collap se (|Forbes et al.lll997n or merg- 
ers (lAshman fc Zepj Il992h . It has also been proposed by 
ICote et al.l ( 19981 ) that red clusters represent the original 
population and blue clusters have been accreted from a 
lower metallicity dwarf galaxy population, in which case 
the low-luminosity ellipticals would be expected to con- 
tain any sma ller population of red, metal-rich clusters 
( Kundu fc Wh itniorc 2001a). In N-body simulations (e.g. 
iMuratov fc Gnedinll2010 ). metallicity bimodality arises from 
the history of galaxy assembly. For instance, early mergers 
are frequently involving relatively low mass protogalaxies, 
which preferentially produce low metallicity blue clusters. 
Late mergers are more infrequent but involve more massive 
galaxies and just a few late massive mergers ca n produce a 
signifi c ant number of red cl usters. Alternative l y. I Yoon et af] 
(|2006l ). lYoon et all (|2011al ). and lYoon et al.l (|2011bl ) sug- 
gested that a broad, single-peaked metallicity distribu- 
tion presumably as a result of continuous chemical enrich- 
ment can produce colour bimodality due to the non-linear 
metallicity-to-colour conversion, without invoking two dis- 
tinct sub-populations. In fact, the study of GC populations 
in low-luminosity ellipticals may be the only way of dis- 
tinguishing between the above models since the predictions 
from all of them are ver y similar for more massive galaxies 
(e.g. lKissler-Pati3l2001^ . It is crucial to test these predic- 
tions for galaxies in both rich clusters and in the field where, 
for example, accretion processes may be by far less efficient. 



The current level of knowledge regarding GC pop- 
ulations in low-l u mino sity ellipticals is quite poor. 
iKissler-Patig et al.l (| 19981 ) concluded on the basis of lim- 
ited ground-based data that the colours of their GC pop- 
ulations were consistent with a unimodal distribution, i.e., 
the colours were very different from those of luminous el- 
lipticals. However, ground-based imaging is severely com- 
promised by the need for large statistical background cor- 
rections. Compounding the problem is the proportionately 
smaller number of clusters associated with lower luminosity 
galaxies. Using compilatio n s of a rchiva l WFPC2 HST data, 
iGebhardt fc Kissler-Pat"i3 l| 19991 ) and iKundu fc Whitmord 
( 2001al ) both found evidence for possible bimodality in a 
small sample of lower luminosity ellipticals. However, the 
limited resolution and depth of the WFPC2 data produced 
discrepancies in the classification of individual galaxies. The 
properties of low-luminosity galaxies in high-density envi- 
ron ments has been add r essed by the ACS Virgo Cluster Sur- 
vey (|Peng et al.ll2006al : IPeng et al.ll2008l ). However, such a 
sample is, by its nature, biased towards environments that 
are atypical of a general field population (Figure [!}. A com- 
plementary field galaxy comparison sample selected from 
low-density environments is required to investigate the role 
of environmental effects on the physical processes that con- 
trol the formation of GC systems and their host galaxies. 



2 SAMPLE SELECTION AND OBSERVATION 

In order to obtain a well-defined sample of early-type galax- 
ies with low-luminosities in low-density environments, a 
complete sample of galaxies was first c omp iled from the 
Nearby Galaxy Catalogue (NBG) (TuUy 1981) with the fol- 
lowing parameters: —5 ^ Hcoda ^ ^fl —18.0 < Mb < 



^ This is a morphological type code. Elliptical and lenticular 
galaxies are included in this range. 
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Figure 1. Distribution of the local density an d the absolute magnitude in B for our sample and the ACSVCS. left panel: The number 
densities are taken from the NBG {Tullv 1988). A clear separation in the environment betwe en our sample (solid histogram) and the 
ACSVCS (dashed histogram) is shown. The distribution of all early-type galaxies in the NBG llTullvl[T98a^ is plotted with a dotted line 
on an arbitrary scale, right panel: The solid histogram represents the distribution of the absolute B magnitude for our sample (Table [TJ, 
while the dashed histogram indicates that of the ACSVCS galaxies. 



-19.^, distance D < 30 Mpc, Gala ctic la t itude |6| > 45 deg, 
and local density po < 1.0 Mpc~^. pTullvl (|l988h determines 
the local density, po, on a three dimensional grid with a 
scale size of 0.5 Mpc after Gaussian smoothing, to estimate 
the contribution of each member of the galaxy population 
brighter than Mb < —16 mag. After a visual inspection 
of the Digitized Sky Survey images, 10 of the brightest (in 
apparent magnitude) targets with no nearby bright stars 
within the ACS field of view were selected. The properties 
of the final sample of early-type galaxies in low-density en- 
vironments are listed in Table [T] in descending order of ab- 
solute luminosity. The local densities of 47 galaxies from 
the ACSVCS sample are also available in the NBG, includ- 
ing all ACSVCS galaxies within the luminosity range of our 
sample. As shown in the left panel of Figure (TJ there is not 
only a clear separation in the local density between our field 
galaxy sample and the ACSVCS sample, but the sample 
in low-density environments is also more representative of 
early-type galaxies in general. The right panel of Figure [1] 
shows the distribution in the absolute magnitude of our sam- 
ple and the ACSVCS sample. 

The observations (Program ID: 10554) were carried out 
with the ACS Wide Field Camera (WFC) on the Hubble 
Space Telescope during Cycle 14 (Oct 05 - Sep 06). The 
WFC is composed of two 4K x 2K chips with a pixel scale 
of 0.05" /pix covering a field of view of 202" x 202" . For each 
galaxy, 2 orbits were used to obtain images in two bands, 

^ The absolute magnitudes of galaxies listed in the NBG scatter 
by ~ 0.5 mag compared to the values from the NASA Extragalac- 
tic Database (NED). We adopt the values from the NED for the 
presents analysis. 



F475W and F850LP, which correspond to Sloan filters g 
and z, respectively. These two filters were also used by the 
ACSVCS, which has the advantage of being able to compare 
results directly without the use of photometric transforma- 
tions. In this work, a four-point line dithering pattern in 
which the telescope pointing shifts 5 x 60 pixels between 
sub-exposures was employed. Such a pattern allows the gap 
between the two chips to be filled and eliminates hot pix- 
els during data processing. Because dithering also effectively 
removes cosmic rays, CR-SPLIT was set to NO. The total ex- 
posure times were maximized by arranging 8 sub-exposures 
in the following way. In the first orbit, four sub-exposures of 
F475W and one sub-exposure of F850LP were allocated. In 
the second orbit, three sub-exposures of F850LP were ob- 
tained. We aim to achieve a S/N ~ 30 at 50% completeness, 
corresponding to estimated exposure times of 23 minutes in 
F475W and 53 minutes in F850LP. The actual total expo- 
sure times depend on the visibility of the target and were 
typically 1300 s and 3000 s for F475W and F850LP, respec- 
tively. The observation logs are listed in Tabled For NGC 
3156, the pointing and orientation were slightly adjusted so 
as to avoid adjacent bright stars. 



3 DATA REDUCTION 
3.1 Basic image reduction 

Each target has four sub-exposure images with different 
dithering positions in each filter. These raw images were re- 
duced using the default ACS pipeline, CALACS, to produce 
bias subtracted and fiat-field corrected images. The default 
processing was found to slightly overestimate the sky back- 
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Table 2. The observation lo 


gs 






Target 


Date 


Filter 


Exp. Time 
(sec) 


Dataset 


NGC 3818 


2006-01-01 


1475W 


1380 


J9CZ04010 






b 850LP 


2987 


J9CZ04020 


NGC 7173 


2006-05-16 


1 475 W 


1375 


J9CZ10U1U 






b 850L1' 


3075 


J9CZ10020 


NGC 1439 


2006-08-21 


b 475 W 


1375 


J9CZUDU1U 






b 850LP 


3023 


J9CZ0d020 


NGC 1426 


2006-08-21 


1 475 W 


1375 


J9CZU9U1U 








3023 


J9CZU9U2U 


NGC 3377 


2006-01-13 


1475W 


1380 


J9CZ08010 






F850LP 


3005 


J9CZ08020 


NGC 4033 


2006-01-04 


F475W 


1380 


J9CZ07010 






F850LP 


3017 


J9CZ07020 


NGC 1172 


2006-08-17 


F475W 


1380 


J9CZ03010 






F850LP 


3005 


J9CZ03020 


NGC 3156 


2005-10-30 


F475W 


1380 


J9CZ02010 






F850LP 


2972 


J9CZ02020 


NGC 3073 


2006-05-14 


F475W 


1428 


J9CZ05010 






F850LP 


3490 


J9CZ05020 


IC 2035 


2006-04-28 


F475W 


1380 


J9CZ01010 






F850LP 


3299 


J9CZ01020 



4 6 8 10 

Radius(pix) 



teertin fc Arnout3ll996h . which models background (galaxy 
light in our case) using interpolation of mean/median values 
in a specified mesh size, thus the galaxy hght in the centre 
may be underestimated. 



Figure 2. Growth curves of bright GCs in the different sample 
galaxies. The dotted lines are growth curves for the ACS PSF. 
The GCs are clearly more extended than a point source. The 
most extended profile is from the closest galaxy (NGC 3377). 



ground because of the large contribution of galaxy light in 
the images. This resulted in negative background levels on 
the outskirts of a galaxy in the final drizzled images. An 
improved sky background was determined by adopting the 
minimum median value from four 200 pix x 200 pix cor- 
ners of each sub-exposure. The four sub-exposure images 
were then combined and geometrically corrected with the 
MultiDrizzle package. The final drizzled image consists of a 
4096 X 4096 pixel science image in units of electrons/s and 
an error map in the second extension that contains all er- 
ror sources such as readout noise, dark current, and photon 
background noise. 

3.2 Galaxy light subtraction 

The smooth galaxy light of each target was fitted with el- 
liptical isophotes using the IRAF/ELLIPSE routine. Fitted 
ellipse parameter tables were then employed in order to build 
a galaxy model image using the IRAF/BMODEL task. The 
model galaxy was subtracted from the original image. It 
was found that subtracting the galaxy continuum light im- 
proves the efficiency of GC detection in the central regions 
of a galaxy. It is because that ellipse fitting does better job 
to model galaxy light in the central regions than SExtractor 



3.3 Globular cluster detection and selection 

SExtractor was run on the galaxy-light subtracted images so 
as to detect GC candidates with DETECT.THRESH > 3a above 
the background. In order to properly account for noise due to 
diffuse galajcy light, the error map from the drizzled images 
was used to create a weight image in SExtractor. Among 
the detected objects, those with ELONGATION > 2 in either 
the F475W or F850LP band and very diffuse objects with 
CLASS_STAR < 0.9 in the F475W band were rejected. Objects 
were matched within a 2 pixel radial separation across the 
two bands. Some spurious detections that were found at the 
edge of the images were manually removed. 

3.4 Globular cluster photometry 

Once the list of GC candidates had been generated, aperture 
photometry was performed on the galaxy-subtracted images 
using the IRAF/PHOT package with 1 to 10 pixel radius 
apertures in one pixel steps. The background was estimated 
locally from a 10 to 20 pixel radius annulus. This locally es- 
timated background has the advantage that the photometry 
of GCs is hardly affected by the global sky background esti- 
mation and the modeling of underlying galaxy light. Because 
GCs are marginally resolved by ACS/WFC at the distance 
of our sample galaxies, the aperture correction from stellar 
PSFs would not be applicable to our GCs. Instead, a rep- 
resentative GC model profile was constructed from several 
moderately bright GCs for each galaxy using the IRAF/PSF 
task. From this process, growth curves were extracted. As 
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Figure 3. The colour-magnitude diagram of GC candidates; the parent galaxies are ordered from most luminous (NGC 3818) to least 
luminous (IC 2035). The dotted box indicates the applied colour and the magnitude cut applied. The number of GCs in this box is 
denoted along with the galaxy name. The background grey scale represents the completeness at a given magnitude and colour. The 
bottom-center panel shows an example of contamination per field, after one out of six objects is randomly selected from 6 compiled blank 
fields. 
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clearly shown in Figure [2l the GCs are not point sources, 
which can be verified by the fact that the light profiles of 
the GCs in the closest galaxy, NGC 3377, are the most ex- 
tended. 

An aperture correction for the GCs in each galaxy was 
determined from the magnitude difference between the 3 
pixel and 10 pixel radius apertures on each growth curve. 
This aperture correction was then applied to the 3 pixel 
radius aperture photometric measurements for all of the 
GC candidates. Such a process can be justified by the 
fact that the sizes of GCs are independent of their lu- 
minosity (I Jordan et al.l |2005| ) . The final corrections from 
the 10 pixel radius to the total magnitudes are primar- 
ily a function of the ACS/WFC PSF and are 0.095 mag 
for F475W and 0.117 mag for F850LP; the zero-points 
of the ABmag scale are 26.068 and 24 .862 for F475W 
and F850LP, respectively (jSirianni et al.l 2005). A correc- 
tion for Galactic extinction was imp lemented based on the 
ISchleeel. Finkbeiner. fc David (jl998f l extinction map values 
(listed in Table [l| and the extinction ratios, Af475w ~ 
3.634^(5 - V) and ^fssolp = 1A85E{B - V), from 
ISirianni et all (|200i). Hereafter, g and z magnitudes refer 
to the F475W and F850LP extinction-corrected total AB- 
mag. Colour-magnitude diagrams for the GC candidates, in 
order of their host galaxy luminosity, are shown in Figure |3l 
with the most luminous galaxy in the top-left panel. 



I I I I I I 




. I I ■ I I ■ ■ 
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Figure 4. The completeness of GC detections in NGC 3377. The 
filled and open circles denote the completeness of the g and z 
bands respectively. The solid and dotted lines are the fitted curves 
from Equation [T] 



Table 3. The parameters of the completeness functions 



Galaxy 



3.5 Globular cluster completeness tests 

To test the effectiveness of the proposed GC detection and 
selection method, a list of 1000 artificial GCs with a uniform 
luminosity function and a spatial distribution in both filters 
was independently generated. These GCs were then added 
onto a galaxy subtracted image that contains all relevant 
noise sources such as the readout noise and poisson errors 
from sky and galaxy light. The artificial GCs were simulated 
using IRAF/MKOBJECTS with the representative GC light 
profile detailed in i]3.4l The GC detection and selection pro- 
cedures described in S I3.3l were then repeated and the initial 
input coordinates of the artificial GCs were matched with 
those of the recovered GCs. For a given magnitude bin, the 
number fractions of recovered GCs were calculated. 

A completeness curve can often be expressed by the 
analytic function. 



F475W 

mo a 



F850LP 
mo a 



f 



mo 



+ a'^{m — mo)'' 



(1) 



where mo is a magnitude when / is 0.5 and a controls how 
quickly / declines (the larger th e value of a, the st eeper 
the transformation from 1 to 0) (|Fleming et al.lll995f ). An 
example of completeness in one of the galaxies is shown in 
Figure U The fitted parameters of Equation [1] are listed in 
Table O The completeness in each galaxy at a given magni- 
tude and colour is displayed as a grey scale in Figure [3l It is 
calculated by multiplying the completeness in g and z over a 
grid of small boxes on the colour-magnitude plane. As seen 
in Figure |4] and Table [S] the completeness levels in g de- 
cline more steeply than in z. This is because the additional 
selection criterion in g, CLASS_STAR, rejects faint GCs more 
quickly due to the increasing uncertainty of CLASS_STAR in 
fainter GCs. 



NGC 3818 


25.94 


4.29 


25.79 


3.06 


NGC 7173 


25.93 


5.24 


25.74 


2.73 


NGC 1439 


26.01 


6.17 


25.84 


3.12 


NGC 1426 


26.05 


5.55 


25.85 


3.44 


NGC 3377 


25.68 


4.81 


25.11 


1.88 


NGC 4033 


25.90 


7.06 


25.81 


3.81 


NGC 1172 


25.86 


3.85 


25.65 


4.03 


NGC 3156 


25.77 


4.88 


25.65 


5.10 


NGC 3073 


25.92 


5.48 


25.79 


3.88 


IC 2035 


26.10 


6.77 


26.02 


4.85 



3.6 Contamination of background galaxies 

Although our selection criteria are quite strict, it is in- 
evitable that some amount of contamination by background 
galaxies will be present (contamination by foreground stars 
is small at these magnitudes with our selection criteria) . To 
estimate this contamination statistically, the ACS archive 
was searched for high Galactic latitude blank-fields (pro- 
posal ID: 9488) that were observed with the same filters 
and were at least as deep as our observations. In total 30 
drizzled images from six blank fields were retrieved. Table 
[4] lists RA, DEC and the number of exposures of each band 
for the chosen blank fields. One field usually consists of 3-6 
sub-exposures that were combined into one image by tak- 
ing median values via IRAF/IMCOMBINE. The aforemen- 
tioned GC detection and selection techniques were applied 
to these images. Since the observation conditions for the 
blank fields are not identical to ours (e.g. exposure time and 
dithering), this difference must be taken into account. It was 
found that some hot pixels were misclassified as GCs in the 
blank fields. These were removed by applying an additional 
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Table 4. Control fields 



RA (J2000) 


DEC (J2000) 


Number of exposures 
for each filter 


14hl8m28.76s 


+24d59m21.95s 


5 


10h34m05.03s 


+58d57m50.60s 


6 


22h22m55.53s 


-72d23ml2.21s 


6 


12h43m30.12s 


+Ild49m20.62s 


6 


12hl0m53.80s 


+39dl4m06.24s 


4 


Ilhl3m33.33s 


+22dl5m41.22s 


3 



classification, FWHM_IMAGE > 1.5. In the bottom-center panel 
in Figure O one out of every six objects in tlie 6 blank fields 
that passed the selection criteria is randomly chosen and 
plotted. Because the blank fields are deeper than ours, the 
completeness of each galaxy was applied so as to obtain the 
colour distribution and luminosity function of the contam- 
ination for each galaxy. These contamination profiles were 
used in the subsequent analysis. The typical contamination 
per field is ~ 14 objects. 



4 RESULTS 

4.1 Colour distributions 

To construct the GC colour distributions, the sample of 
GC candidates was refined by applying the colour cut, 
0.6 < {g — z)o < 1.7, and the magnitude cut, Mg > —12.0, 
and g < 25.5, where the completeness is ~ 80%. These 
colour and magnitude selections are indicated by the dotted 
boxes in Figure [31 In Figure [S] the raw distribution is drawn 
with a black histogram. The red histogram represents the ex- 
pected contamination level mentioned in the above section. 
For each colour bin, the number of contaminating objects 
was removed from the raw distribution; the green histogram 
represents the contamination-corrected colour distribution. 
The panels in Figure [S] are placed in the same order as in 
Figure [S] 

In order to statistically test the significance of any 
colour bimodality, the Kaye 's mixture model (KMM) test 
jAshman. Bird, fc Zeplll994l ') was employed. The KMM test 
uses likelihood ratio test statistics to estimate the probabil- 
ity (P-value) that two distinct Gaussians with the same dis- 
persion are a better fit to the observational data than a sin- 
gle Gaussian. After running the KMM test for 100 bootstrap 
resamples of the contamination corrected colours, the peaks 
of the colour histograms, the a of the sub-populations, and 
the red population fraction were estimated. Also estimated 
were the P-value and the fraction of P-values that are less 
than 0.05 (Table[5la)); the representative values are median 
values of the 100 KMM outputs and the uncertainties are 
half a width within which 68% of the data are contained rel- 
ative to the median. Since the KMM outputs (especially the 
P-values) have a somewhat skewed distribution, the median 
is a more robust estimation than the mean. The statistical 
interpretation of the P-value is that a distribution with a P- 
value of 0.05 favors bimodality rather than unimodality with 
a 95% confidence. In Figure [S] the results of the KMM tests 
for each subgroup are plotted with dashed lines, and the sum 
of the two groups is plotted with solid lines. The significance 



of colour bimodality is decided on the basis of (P) and the 
fraction of P < 0.05. If (P) < 0.05 and /(p<o.o5) > 0.90, a 
distribution is interpreted as a strong bimodal distribution 
and for 0.05 < (P) < 0.10 or 0.50 < /(p<o.o5) < 0.90, a 
distribution is likely bimodal; otherwise it is unimodal. De- 
spite the small P-value for IC 2035, this galaxy is regarded as 
unimodal because of the insignificant number of red clusters. 
The significance of bimodality is listed as S(trong), L(ikely), 
or U(nimodal) in Table [S] 

The ability of the KMM test to detect bimodality de- 
pends on the sample numb er and the norma lized separation 
between the two peaks (see I Ashman. Bird, fc Zcpf 1994). A 
larger sample size and a larger separation between two sub- 
components will increase the chance of detecting bimodality. 
Furthermore, the extended tails of the distribution can also 
affect the results. Therefore, an alternative GC sample with 
a narrower colour range of 0.7 < {g — z)o < 1.6 was tested. 
The corresponding results are listed in Table [S{b). Clipping 
the tails of the distributions was found to decrease the P- 
value. In other words, without the extended tails, the KMM 
test is more likely to detect bimodality. This is not surpris- 
ing because fitting extended tails results in a larger sigma for 
each subgroup, thus making it harder to separate two sub- 
groups. However, in this experiment, the values of most of 
the fitted parameters do not vary significantly and remain 
within the uncertainty that was initially estimated. There 
is one case where a fitted parameter, the red peak colour, 
has been changed significantly after clipping the tails. This 
scenario occurs in NGC 3156, where one isolated and very 
red GC (see middle-bottom panel in Figure [S]) makes the 
red peak significantly redder compared to the value when 
it is removed. This is potentially a background galaxy con- 
taminant because, unlike other red GCs, it lies at a large 
projected distance from the galaxy center. For further anal- 
ysis, the values from Table (5^ a) were adopted in this study 
except for NGC 3156, where the value from Table \^h) is 
used. The KMM results are overplotted in Figure [5] 

It should be noted that colour bimodality in GCs is 
more common in the more luminous galaxies. The red pop- 
ulation becomes weaker and moves to bluer colours with 
decreasing galaxy luminosity. Eventually, only a blue peak 
appears in the faintest galaxies. The richness of the GCs 
also decreases with the host galaxy luminosity. None of our 
GC colour distributions have a single broad peak located 
between the normal blue and red peaks. This is in contrast 
to the results of so me studies on early-type galaxies in other 
environments (e.g. Larsen et al.l 1200 ll : iKundu fc Whitmord 
l2001al : iPeng et al.ll2006al '). The behavior of the GC colour 
with host galaxy luminosity will be discussed in more detail 
below. 



4.2 Luminosity functions 

Luminosity functions in the g and z bands were constructed 
using GCs in the colour range of 0.6 < (<; — z)o < 1.7 (black 
histograms in Figure [6] and Figure [T]). The contamination 
(red histograms) estimated previously was subtracted from 
the raw distribution and incompleteness corrections were ap- 
plied. The corrected luminosity functions (green histograms) 
were fitted with a Gaussian function up to the 50% com- 
pleteness limit. The fitted parameters (peak magnitude and 
standard deviation) and their uncertainties are listed in Ta- 
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Figure 5. Colour distributions of GC candidates. The black histograms are colour distributions of the raw data. The red histograms are 
contamination estimates. The green histograms are after correcting for contamination. The KMM results for sub-populations are drawn 
with dashed lines, and the solid lines are the sum of the two sub-populations. 
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Figure 6. Luminosity functions of GCs in the g band. The blacls histograms are luminosity functions before correction. The red histograms 
are contamination estimates, and the green ones are luminosity functions corrected for background contamination and completeness. 
The solid curves are the best fit Gaus sian luminosity funct ions, and the arrows indicate the expected turn-over position for each galaxy 



1 ne sona curves are tne best nt Lraus sian mmmosity lunct i 
based on the GCLF parameters from lljordan et al.ll2007bl) 




Figure 7. Luminosity functions of GCs in the z band. The histogram eolours and fitted curves have the same meaning as for the g band 
luminosity functions in Figure [6] 
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Table 5. The properties of globular cluster systems: colour distributions 



Host galaxy 


Ngc 


Blue Peak 


Red Peak 


suh 


fred 


(P) 


/(P < 0.05) 


Bimod. 


{(9 


- z)o) 


(a) 














0.6 < {g 


-z)o < 1.7 














NGC 3818 


240 


0.95 


± 


0.01 


1.33 ± 0.02 


0.12 


0.46 


0.00 ± 0.00 


1.00 


S 


1 


12 


± 


0.01 


NGC 7173 


208 


0.96 


± 


0.01 


1.38 ± 0.02 


0.12 


0.40 


0.00 ± 0.00 


1.00 


S 


1 


13 


± 


0.02 


NGC 1439 


139 


0.98 


± 


0.03 


1.32 ± 0.04 


0.15 


0.35 


0.02 ± 0.09 


0.63 


L 


1 


10 


± 


0.02 


NGC 1426 


159 


0.93 


± 


0.02 


1.27 ± 0.02 


0.12 


0.37 


0.00 ± 0.00 


0.99 


S 


1 


06 


± 


0.02 


NGC 3377 


173 


0.94 


± 


0.02 


1.30 ± 0.02 


0.13 


0.52 


0.00 ± 0.01 


0.96 


s 


1 


13 


± 


0.02 


NGC 4033 


111 


0.94 


± 


0.03 


1.29 ± 0.03 


0.14 


0.38 


0.03 ± 0.16 


0.56 


L 


1 


07 


± 


0.02 


NGC 1172 


265 


0.96 


± 


0.01 


1.24 ± 0.02 


0.09 


0.32 


0.00 ± 0.00 


1.00 


S 


1 


05 


± 


0.01 


NGC 3156 


74 


0.95 


± 


0.05 


1.29 ± 0.20 


0.13 


0.16 


0.03 ± 0.09 


0.57 


L 


1 


01 


± 


0.02 


NGC 3073 


52 


0.98 


± 


0.03 




0.13 


0.06 


0.28 ± 0.42 


0.16 


U 


1 


02 


± 


0.02 


IC 2035 


35 


0.97 


± 


0.03 




0.12 


0.06 


0.00 ± 0.01 


0.83 


U 


1 


00 


± 


0.03 


(b) 














0.7 < {g 


-z)o < 1.6 














NGC 3818 


236 


0.95 


± 


0.01 


1.33 ± 0.01 


0.11 


0.46 


0.00 ± 0.00 


1.00 


s 


1 


12 


± 


0.01 


NGC 7173 


201 


0.97 


± 


0.01 


1.37 ± 0.02 


0.11 


0.40 


0.00 ± 0.00 


1.00 


s 


1 


14 


± 


0.02 


NGC 1439 


134 


0.98 


± 


0.02 


1.30 ± 0.04 


0.13 


0.39 


0.01 ± 0.02 


0.84 


L 


1 


10 


± 


0.02 


NGC 1426 


157 


0.94 


± 


0.02 


1.28 ± 0.02 


0.11 


0.38 


0.00 ± 0.00 


1.00 


s 


1 


07 


± 


0.02 


NGC 3377 


165 


0.95 


± 


0.01 


1.30 ± 0.01 


0.11 


0.53 


0.00 ± 0.00 


1.00 


s 


1 


13 


± 


0.02 


NGC 4033 


105 


0.95 


± 


0.01 


1.30 ± 0.02 


0.11 


0.40 


0.00 ± 0.00 


0.98 


L 


1 


09 


± 


0.02 


NGC 1172 


264 


0.96 


± 


0.01 


1.24 ± 0.01 


0.09 


0.33 


0.00 ± 0.00 


1.00 


s 


1 


05 


± 


0.01 


NGC 3156 


73 


0.92 


± 


0.03 


1.20 ± 0.05 


0.11 


0.23 


0.09 ± 0.15 


0.40 


L 


1 


00 


± 


0.02 


NGC 3073 


50 


0.99 


± 


0.03 




0.11 


0.16 


0.08 ± 0.27 


0.42 


U 


1 


03 


± 


0.02 


IC 2035 


34 


0.98 


± 


0.02 




0.11 


0.06 


0.00 ± 0.00 


0.88 


u 


1 


01 


± 


0.03 



Table 6. Gaussian fitted parameters of luminosity functions and specific frequencies 



Host galaxy 


mto, 


9 








mto,z 




Sn 


NGC 3818 


25.09 




0.13 


1.05 


± 


0.09 


24.02 ± 0.17 


1.14 ± 0.11 


1.36 ± 0.12 


NGC 7173 


25.28 


± 


0.15 


1.04 


± 


0.11 


24.36 ± 0.24 


1.12 ± 0.14 


1.63 ± 0.06 


NGC 1439 


24.66 


± 


0.08 


0.77 


± 


0.07 


23.51 ± 0.10 


0.84 ± 0.08 


0.91 ± 0.10 


NGC 1426 


24.73 




0.12 


1.02 


± 


0.09 


23.77 ± 0.14 


1.11 ± 0.11 


1.67 ± 0.21 


NGC 3377 


23.08 


± 


0.12 


1.25 


± 


0.11 


21.75 ± 0.11 


1.14 ± 0.09 


2.35 ± 0.49 


NGC 4033 


24.88 


± 


0.24 


1.22 


± 


0.17 


24.01 ± 0.33 


1.26 ± 0.20 


1.91 ± 0.19 


NGC 1172 


24.18 




0.09 


1.06 




0.07 


23.19 ± 0.07 


0.96 ± 0.06 


9.18 ± 4.41 


NGC 3156 


















0.62 ± 0.02 


NGC 3073 


25.17 


± 


0.67 


1.53 


± 


0.44 






3.48 ± 1.76 


IC 2035 


24.60 


± 


0.19 


0.79 


± 


0.17 






0.91± 0.15 



ble [6l the missing data correspond to parameters with a 
large uncertainty caused by the small number of GC sam- 
ples and an excess of faint GCs. For the luminous galaxies, 
Gaussian functions aptly represent the luminosity functions 
of the GC systems. The arrows in Figures [5] and [7] indicate 
the expected turn-over magnitudes based on a peak abso- 
lute magnitude of the GC luminosit y function (GCLF) o f 
= —7.2 and [i,^ = —8.4 as given by Jordan et al.l (|2007bl ). 

As can be seen in Figure [S] there is an excess of faint 
GCs around g ~ 25.5 mag in the GCLF of NGC 3377. 
These faint GCs deviate significantly from the fitted Gaus- 
sian function. This deviation is also found at z ~ 24.5 mag 
in the z-band GCLF shown in Figure [T] The properties of 
these objects are discussed in more detail below. 

In Figure [S] the spatial distribution, size, and colour of 
NGC 3377 GCs fainter than g = 25 mag are plotted along 
with those brighter than g = 25 mag. In the top-left panel. 



the faint GCs (large closed circles) appear to be uniformly 
distributed, unlike the bright GCs (small open circles) which 
are concentrated at the galaxy center. Shown in the top-right 
panel of Figure 8 are histograms of the FWHM returned by 
Sextractor for the faint and bright GCs. The size distribution 
of the faint GCs is also quite different from that of the bright 
GCs in the sense that the sizes are either smaller or larger 
than the medium size for all clusters. The colour distribution 
of the faint GCs appears to be rather flat, as shown in the 
bottom-left panel. Any dependence of the colour on the size 
of the faint GCs was examined by separating the GCs at 
FWHAI = 2.8 pixels. In the bottom-right panel of Figure 
8, the thin histogram represents the colour distribution of 
faint GCs with FWHM > 2.8, while the thick histogram 
represents GCs with FWHAd < 2.8. No obvious correlation 
between size and colour was found in the plot. 

We can speculate on the possible origin of the faint 
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Figure 8. Properties of faint GCs near g ~ 25.5 mag in NGC 3377. The top-left panel shows the spatial distribution of bright GCs (g < 25 
mag, small open circles) and faint GCs {g > 25 mag, large closed circles). The top-right panel shows the histograms of FWHM of the 
bright (thin histogram) and faint (thick histogram) GCs. The bottom-left panel shows colour distributions of the bright (thin histogram) 
and faint (thick histogram) GCs. The bottom-right panel shows colour distributions of the faint GCs with larger FWHM > 2.8 (thin 
histogram) and smaller FWHM < 2.8 (thick histogram). 



objects tha t devia te from the normal Gaussian GCLF. 
IPeng et all (l2006b') found a class of diffuse star clusters 
(DSCs) in the ACSVCS and investigated the nature of these 
objects, which are characterized by their low- luminosity, 
broad distribution of sizes, low surface brightness, and red- 
der colour when compared to normal metal-rich GCs. The 
spatial distribution of the objects was also closely associated 
with the host galaxy light. However, not all of the aforemen- 
tioned characteristics are found in our faint objects, which 
possess a wide distribution of sizes, but a significant number 
of blue GCs and a random spatial distribution. 

Another possibility is that the faint objects may be 
field stars (bright giants) that belong to the parent galaxy. 
By comparing isochrones of various ages and metallicities 
l|Girardil [20061 ). the locus of giant stars {My < —5) of 
a very young age (~10 Myr) was found to overlap with 
that of the faint objects on the colour-magnitude diagram 
(Figure However, it is quite unlikely that such young 
stars exist in significant numbers in old elliptical galaxies 
like NGC 3377. Furthermore, horizontal branch stars are 



not sufficien t ly lum inous to account for the faint objects. 
[Harris et al.l (|2007l ) found no evidence of young stellar pop- 
ulations (< 3 Gyr) in the halo star CM diagram of NGC 
3377 from deep HST/ACS photometry, although their field 
of view barely overlaps that used in this study. The number 
density of foreground Milky Way stars in the direction of 
NGC 3377 is also very low {Nstar < 10) within the ACS 
field (|Robin et al.ll2003l ) and even lower in the magnitude 
range of the faint GCs. Visual inspection of the faint objects 
reveals that some are more likely to be misclassified back- 
ground galaxies. Assuming the faint objects are not GCs 
belonging to NGC3377 shifts the mean colour of the GC 
system blueward by {A{g — z)q 0.015) compared to the 
value including all the objects. Note that the mean colours 
in Table [5] have been estimated after a faint magnitude cut 
{g = 25.5 mag) so that any colour change by faint objects is 
negligible. 
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Figure 9. A plot of GC colours against host galaxy luminosity. 
The diamond and circle symbols represent the blue and red peaks 
respectively. Strong bimodality is indicated by a filled symbol, and 
marginal or zero bimodality by an open symbol. The mean colours 
of all GCs within a galaxy are plotted with an open triangle. The 
two solid lines are the least squares best fits to the blue and red 
peaks. The dotted line is fitted to the mean colour. 




Figure 10. Diff'erences in g and z magnitudes and g — z colour for 
GCs in 24 ACSVCS g alaxies betwe e n our aperture photometry 
and KING PHOT from ljordan et al.l ||2009|) . The measurement of 
our aperture photometry is on the x-axes and offsets between the 
two measurements are on the y-axes. From top to bottom the 
offsets plotted are: Ag = gap — gxing , = Zap — ZKing and 

A(g - z)o = (g - z)o,ap - (g - z)o,King- 



4.3 Host galaxy properties vs. globular cluster 
colours 

It is well known that there is a strong correlation between 
mean colours /red p eaks and host galaxy lumin osities (e.g. 
iLarsen et al] I2001I : iKundu fc Whitmorel l2001ah . However, 
few studies had been performed to support a correlation 
between the blue peak and the host galaxy luminosity un- 
til the recent ACSVCS data found such a relati onship with 
their l arge dynamic range of galaxy luminosity (|Peng et al.l 
l2006al ). 

In Figure|9l the peak colours of the sub-populations and 
the mean GC colours are plotted against the host galaxy 
luminosity for our sample. Each relationship is fitted to a 
straight line using weighted chi-squre minimization. It is 
clear that strong dependencies of the red peak and the mean 
GC colour on the host galaxy luminosity do exist in our field 
galaxies with almost identical slopes. On the other hand, the 
blue peak is almost independent of the galaxy luminosity. 
In fact, the fitted slope exhibits a weak anti-correlation, the 
statistical significance of which is within la of a zero slope. 
From the attained results alone, it is uncertain whether the 
absence of a blue peak gradient is due to the narrow dy- 
namic range of the luminosity in our field galaxy sample or 
an intrinsic effect across a wide range of luminosity. The ob- 
tained results clearly show that the observed trends in the 
GC colours with host galaxy luminosity are a consequence 
of the gradual change in the fraction and colour of the red 
GCs and not a systematic change in the colours of all GCs 
within a galaxy. 



4.4 Comparison with the ACS Virgo Cluster 
Survey 

The ACS Virgo Cluster Survey (ACSVCS) was a large pro- 
gram used to image 100 early- type galaxies in the Virgo clus- 
ter w ith the ACS in two filters {g and 2 bands) ijCote et al.l 
'2OO4I). The main scientific goals were to study the properties 



of GC systems in these galaxies (e.g., PJC0 6. iJordan et al 



^06), analyze their central structures (jFerrarese et al 



2006h . an d obtain accura te surfa ce brightness fiuctuat ion 
distances (|Mei et al.ll2007l : see also iBlakeslee et al.llioogi '). 

The observation conditions and data reduction proce- 
dures of the ACSVCS are somewhat different from those 
adopted in this study. The ACSVCS used 100 orbits of the 
HST, allocating only one orbit for each galaxy. Each galaxy 
of the ACSVCS had an exposure time of 750 s in F475W and 
1210 s in the F850LP filter, while our exposure times are al- 
most twice as long (typically 1300 s in F457W and 3000 s in 
F850LP). The data reducti on procedures of th e ACSVCS 
are described in detail by Ijordan et al.l l|2004h . In brief, 
the exposures were split into two or three sub-exposures, 
without repositioning, so as to remove cosmic rays. These 
sub-exposures were then reduced and combined using the 
standard ACS pipeline. Model galaxies we re created using 
the ELLIPROF program ijTonrv et al.ll 19971 ) and subtracted 
from the original images. The KINGPHOT program devel- 
oped by Jordan et al. (2005) was used for photometry of 
the GCs and to measure their sizes. KINGPHOT fits King 
models convolved with a given point spread function in each 
filter to individual GC candidates and finds the best fit pa- 
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Figure 11. Comparison of mean colours for GCs in our 
galaxy sample with those from the ACS Virgo Cluster Survey 
(ACSVCS). Our results are plotted with black filled triangles, 
while ACSVCS results from PJC06 are plotted with grey open 
triangles. The mean colours of our field galaxy sample appear to 
be situated at the bluer end of the colour range for a given galaxy 
luminosity. 




Figure 12. Comparison of colour peaks with the ACSVCS. Sym- 
bols are identical with Figure (9] Black symbols are from our 
results, while grey symbols are from the ACSVCS. The dark 
sold lines are a linear fit to each sub-population from our re- 
sults, whereas the grey solid lines are linear fits by PJC06 to the 
ACSVCS. 
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Figure 13. Comparison of fraction of red GCs with the ACSVCS. 
The black filled and open circles represent strong and likely bi- 
modality from this work respectively, while open circles are from 
the ACSVCS. No corrections have been made for any aperture 
bias effects between the two samples (see text for details). 



rameters of the King model. Total magnitudes for the GCs 
were calculated by integrating these best-fit convolved King 
models. All GCs were selected based on their magnitude 
{g 19.1 ox z ^ 18.0) and mean elongation ((e) ^ 2, where 
(e) = a/6) in the two filters. The main results regarding the 
colour distributions of the ACSVCS GC systems were pub- 
lished by PJC06 and are compared with the results of this 
study below. 

Before directly comparing results from the ACSVCS 
with ours, it is necessary to check consistency of photom- 
etry, because the applied photometry methods are differ- 
ent; aperture photometry has been used in this work, while 
the ACSVCS used fitted King models. We applied our pho- 
tometry procedure to 24 galaxies of the ACSVCS within 
our galaxy luminosity range and directly compared mag- 
nitudes and colours of GCs we mea s ured with those in 
the GC catalogue from [Jordan et all (|2009l ). Their cata- 
logue lists two kinds of photometry values, PSF convolved 
King model and aperture photometry. Since PJC06 used 
the former for their analysis, it was taken here for com- 
parison. Figure [TU] shows differences of g and z magni- 
tudes (Ag = gap - gKing,Az = Zap — ZKing) and ColoUr 

(A(g - z)q = {g - z)o^ap - {g - z)o,King) for aU the GCs 
matched in 24 ACSVCS galaxies. It is found that mean 
magnitude offsets in g and z band are (Ap) ~ 0.025 and 
{Az) ~ 0.024 respectively, revealing that our photometry 
appears to slightly overestimate total flux of GCs. However, 
no offsets in mean colour between the two methods are seen, 
{A{g — z)o) < 0.01. This test verifies that there is no sys- 
tematic colour offset between the two methods so that direct 
comparison of GC colours between PCJOG's results and our 
results can be justified. 
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In Figure [TT] the mean colours for the 10 GC systems in 
our sample are compared with the results from PJC06. At a 
given host galaxy luminosity, our mean colours are generally 
located at the bluer end of the mean colour distribution of 
the ACSVCS galaxies. However, unlike the other GC sys- 
tems in our sample, the mean colour of the NGC 3377 GC 
system lies in the middle of the range of ACSVCS mean 
colours. Because the distribution of red GCs is more cen- 
trally concentrated and NGC 3377 is the nearest galaxy in 
our sample, it is possible that the GCs detected within our 
ACS field are biased toward red GCs. To test whether this 
selection bias is large enough to shift the mean colour, the 
GCs in other galaxies were resampled with a smaller field of 
view and their mean colours were recalculated. For this test 
NGC 4033 and NGC 1172, which have a similar absolute 
magnitude but are as twice as distant as NGC 3377, were 
chosen. The mean colour of the GCs in these two galax- 
ies within a 50" radius from their galaxy centers is ~ 0.03 
mag redder than our original estimations for the full spatial 
coverage. Noting that the full field of view of the ACS is 
202" X 202" so that a circle with a 50" radius covers almost 
half of the original aperture, the intrinsic mean colour of the 
GCs in NGC 3377 could be about ~ 0.03 mag bluer than 
our original estimation. However, this is a relatively small 
correction when compared to the wide range of mean colours 
observed in the ACSVCS. Note that the error bar of NGC 
3377 plotted in Figure [TT] is also ±0.02 mag. 

PJC06 analyzed their colour distributions with two 
different methods. One method separates the two sub- 
populations in each host galaxy and tests for bimodality us- 
ing the KMM routine. The other method coadds the colour 
distributions of GCs in bins of the host galaxy luminosity 
and separates the two sub-populations by nonparametric de- 
composition. The results from the first method of PJC06 
were compared with our results because our colour distri- 
butions were also analyzed with the KMM method and the 
number of our sample galaxies is not large enough or wide 
enough to bin by galaxy luminosity. A comparison of our 
results and those of PJC06 with respect to the colour peaks 
of the two sub-populations is shown in Figure 1121 A com- 
parison of the fractions of red population GCs detected was 
also plotted; the results are shown in Figure [131 In Figure [T3l 
only galaxies with a bimodal colour distribution are plotted. 

Overall trends for the colour distributions with respect 
to the host galaxy luminosity in both samples appear to be 
similar. As evident in Figures [T^ and [T51 the probability of 
bimodality decreases with decreasing host galaxy luminosity 
as the red population becomes weaker. There is a strong cor- 
relation between the peak of the red sub-population and the 
host galaxy luminosity. Moreover, the fraction of red GCs 
with respect to the entire population appears to increase 
with host galaxy luminosity. 

Some differences do exist between the results of this 
study and those of PJC06. In Figure 1121 our slope for the 
red peak against galaxy absolute magnitude appears to be 
steeper than that of the ACSVCS. In other words, the bi- 
modal colour distributions in our sample disappear more 
quickly as galaxies become fainter. Furthermore, we found 
no correlation between the blue peaks and the galaxy abso- 
lute magnitude, while PJC06 found a weak correlation (al- 
though it was not as steep as for the red peaks). Our fraction 
of red GCs was also found to be slightly lower than that of 



the ACSVCS at a given galaxy luminosity. The mean colours 
of the GCs in our sample are slightly bluer than those of the 
ACSVCS and the red and blue peak positions in our sample 
are almost identical to those in the ACSVCS. This suggests 
that the difference in the mean colours is due to differences 
in the relative fraction of the red population and not an 
overall shifting of colours. In the middle bottom panel of 
Figure [14] it appears as if the ACSVCS galaxies have more 
red objects. To test whether this is real or not, a K-S test 
was performed, but its results give a low significance level 
to any difference. 

One might argue that because the majority of our sam- 
ple galaxies are more distant than the Virgo cluster galaxies, 
and red GCs are more concentrated on the galaxy center, our 
detected GC samples would be less biased toward red GCs. 
Consequently, our red GC fraction could be systematically 
lower than that in the ACSVCS. Indeed NGC 3377, the near- 
est galaxy in our sample and closer than the Virgo cluster, 
has the highest red GC fraction. To test this hypothesis, we 
resampled GCs in NGC 3818, which is as twice distant as 
Virgo, within a circle with 50"radius that covers half of the 
ACS field of view. Therefore, the actual physical overage of 
this circle is similar to that of the ACSVCS. The new esti- 
mation of the red GC fraction within this circle is slightly 
higher by ~0.02 mag than our initial estimate. In the case 
of NGC 7173, the original image contains 3 galaxies (NGC 
7173, NGC 7174 and NGC 7176) and the actual coverage 
of NGC 7173 used for GC detection is only half of the full 
image. Considering NGC 7173 is as twice distant as Virgo, 
both physical coverages are similar, but the red fraction of 
GCs in NGC 7173 is stiU lower than that of the ACSVCS 
at the same galaxy luminosity. We therefore conclude that, 
although the ACSVCS samples are possibly biased to red 
GCs due to smaller physical coverage, the derived red GC 
fraction does not shift significantly and we still find a small 
systematic difference with the ACSVCS which may be re- 
lated to the effects of local environment. 

The shapes of the colour distributions in the ACSVCS 
are more varied than those in our sample. For instance, 
PJC06 found one single broad peak for the colour distri- 
butions of VCC 1664 and VCC 1619. They also found ex- 
amples in which the red population completely dominates; 
red GC fractions of 0.84 in VCC1146 and NGC 4458 were 
measured. These types of "abnormal" colour distributions 
are not found in our smaller field sample. 

To summarize the comparison of the properties of the 
GC systems in our sample with those of the ACSVCS, plots 
of the GC properties against host galaxy absolute magni- 
tude (identical to Figures 1111 1121 and I13|) are shown in the 
first column of Figure 1141 The second and third columns in 
this figure are comparisons of the GC properties against the 
local galaxy number density (environment) and morpholog- 
ical type, respectively. In these two columns, only results 
from the ACSVCS within the range of the galaxy luminos- 
ity of our sample ( — 18.4 Mg ^ —20.4) are plotted so as 
to minimize any galaxy luminosity dependence. In the third 
column, the GC systems with the same host galaxy morpho- 
logical type are randomly positioned within the morpholog- 
ical type bins. As can be seen in the second column, our 
sample is clearly distinct from the ACSVCS in terms of the 
galaxy environment. Because the GC system properties of 
both the ACSVCS and our sample are widely spread how- 
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Figure 14. Comparison of tlie global properties of GC systems in our sample of low-luminosity field E/SOs with those of the ACSVCS. 
The first, second, and third columns show plots versus the absolute magnitude, the local density of galaxies from TuUy (1988), and the 
morphological type of the host galaxy, respectively. The rows show the colour of the blue and red peaks, mean colour, and the fraction 
of red GCs. The coloured symbols represent our data, while the grey symbols are from the ACSVCS. For the second and third columns, 
only galaxies from the ACSVCS within the luminosity range of our sample are plotted. For the first and third rows, closed circles show 
the GCs in which colour bimodality is strongly detected using the KMM test, while open symbols show GCs with weak colour bimodality. 
In the top-left panel, black lines are the best linear fit to our data and the grey lines are the best-fit to the ACSVCS data. 
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Figure 15. Gaussian fit parameters of luminosity functions against host galaxy luminosity. Turn-over magnitudes (left) and dispersions 
(right) in the g and z bands against h ost galaxy luminosity. The open circles represent the Gaussian fits to galaxy-binned GCLFs from 
the ACSVCS by I Jordan et al! I ll2007d) . whereas the results from this work are plotted with filled circles. 



ever, distinct differences with respect to the galaxy number 
density are not obvious in these plots. In the plots of the 
GC properties against morphological type (third column of 
Figure I14|l no obvious trends are observed with morphologi- 
cal type for the ACSVCS. However, for our sample, the GC 
systems in SBO/SO type galaxies appear to have little/no 
bimodal colour distributions. In other words, the SBO/SO 
galajcies have smaller red sub-populations when compared 
to early- type galaxies. 

Our GCLFs have also been compared with those of 
the ACSVCS and the results are shown in Figure 1151 
[Jordan et aD ()2007b fl fitted the GCLFs of 89 galaxies with 
both a Gaussian function and a Schechter function. They 
found that the GCLF dispersions were correlated with the 
galaxy luminosity, while the GCLF turn-ov er magnitudes 
were r ather const ant (in both bands). Later IVilleeas et al.l 
l|2010l ) reanalvsed [Jordan et al.l l|2007bl )'s results for the For- 



nax cluster, finding that the turn-over magnitude of GCLF 
gets fainter with decreasing galaxy luminosity. Using dis- 
tance moduli from Table [T] and apparent turn-over magni- 
tudes from Table [6l we directly compared our results from 
i|4.2l with the Gaussian fit ting parameters of t he galaxy- 
binned GCLFs reported by|jord|nelaD (|2007bl ). In Figure 
15\ our result s are ove rplotted on the fitting parameters from 
Jordan et all l|2007bl ) and appear to be more widely spread 
than those from the ACSVCS. However, one must keep in 
mind that the ACSVCS results come from galaxy-binned 
GCLFs, which are less noisy than our individual GCLF fit- 
ting parameters. Using only our GCLF dispersions, it is hard 
to say whether or not there is a trend with galaxy luminosity 
because of the narrower range of galaxy luminosity in our 
sample. Our GCLF dispersions agree well with the ACSVCS, 
at least within the range of galaxy luminosities spanned by 
our sample. 
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Figure 16. Evolution of g — z colour for various metallicities 
and ages. The grey d ashed tracks ar e taken from simple stellar 
population models bv iGirardil l l2006h . while the black ones are 
from YEPS. The metallicity noted at the end of each track is in 
units of [Fe/H]. The two stripes represent the colour ranges of 
red and blue peaks of GC systems in our sample galaxies. For old 
stellar populations (> 8 Gyr), {g — z)o colour is mainly governed 
by metallicity rather than age. 



Figure 17. Metallcity and age differences of GC systems between 
high density (cluster) and low density (field) environments. The 
black solid, dotted, dashed and dot-dashed lines represent 10, 11, 
12 and 13Gyr of cluster GCs from the YEPS models respectively. 
The grey solid, dotted, dashed and dot-dashed lines are 10, 11.2, 
12.6 and 14.1 Gyr models from GIR06 respectively. These model 
lines were derived from our results that the mean colours of GC 
systems in cluster galaxies are 0.05 mag redder than in field galax- 
ies. 



5 DISCUSSION 

5.1 Effect of environment on galaxy formation 

From the analysis presented above, it appears that the 
mean colours of GC systems in field environments are 
slightly bluer than those in cluster environments at a 
given host galaxy luminosity. The simplest interpretation 
of this finding would be that these GC systems are ei- 
ther less metal-rich or younger than their counterparts in 
rich clusters (or both metal-poor and young). Simple stellar 
population models from Giraradi (2006, hereafter GIR06) 
and the Yonsei Evolutionary Population Synthesis (YEPS, 
[http://web . yonsei. ac.kr/cosmic/data/YEPS. htm ; C. Chung 
et al., in prep.) model were adopted in order to estimate the 
magnitude of these effects and compare them with our mean 
colours. The prediction of {g — z)o colour in various age and 
metallicity ranges based on the GIR06 and YEPS models is 
shown in Figure 1161 The mean metallicity of the blue and 
red peaks shows good agreement with predictio ns based on 
the empirical transformat ion (g ~ z) to [Fe/H] (|Peng et al.l 
l2006al : ISinnott et ai]|2010l ) . 

With only one colour {g — z), it is almost impossible 
to break the well-known age-metallicity degeneracy. Thus, 
an attempt was made to find any combination of metal- 
licity and age differences that reproduce the mean colour 
offset between the high- and low-density GC systems. We 
first assumed an age range from 10 Gyr to 14 Gyr for the 



cluster environment GC systems and a minimum age of 8 
Gyr for the field environment GC systems. We then took the 
mean colour of the cluster and field GC systems from Fig- 
ure[TT]and derived the metallicity difference for a certain age 
combination from the YEPS and GIR06 models. The metal- 
licity and age differences for cluster GC systems of various 
ages is shown in Figure [TT] As an example, if the mean age 
of the cluster GC systems is 10 Gyr (the black solid line 
in Figure I17p and that of the field GC systems is 8 Gyr, 
then there seems to be no metallicity difference between the 
two GC systems according to the YEPS model. However, if 
the cluster and field GC systems are coeval, the cluster GC 
systems are, on average, more metal-rich by ~ 0.10 — 0.15 
dex than the field GC systems. The YEPS model predicts 
a slightly lower metallicity difference than the GIR06 model 
when the two samples are coeval. This metallicity offset 
between the two samples disappears if the GC systems in 
the field galaxies are younger than those in the clusters. 
The age differences with no metallicity offset depend on the 
assumed ages of the GC systems in high-density environ- 
ments and on the adopted simple stellar population model. 
However, the colour offset can be successfully explained by 
at least a ~ 2 Gyr age offset when there is no metallicity 
difference. Therefore, we tentatively conclude that GC sys- 
tems in low-density environments are either less metal-rich, 
A[Fe/H] ~ 0.10 - 0.15, or younger, Aa^e > 2 Gyr, than 
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those in high-density regions. Of course, any combination of 
age and metallicity within the above ranges can be a possible 
solution, as seen in Figure [T7l 

Metallicity difference between the low- and high-density 
GC systems has been also deri ved from a ful l y em piri- 
cal colour-metallicity relation by iBlakeslee et al.l (|2010l ) , in 
which a quartic polynomial is fit to compiled data from 
PJC06. Based on the relation, the colour offset is found to 
correspond to /S.[Fe/H] ^ 0.12. Although the empirical re- 
lation might contains age effects, the metallicity difference 
is well consistent with that from the SSP models. 

Even though we managed to detect the colour off- 
set of GC systems in low- and high-density environments, 
and quantify the age and metallicity difference using the 
simple stellar population models and the empirical colour- 
metallicity relation, the detected colour offset {A{g — z)o ~ 
0.05) is not comparable with the colour range due to host 
galaxy luminosity. For example, within our sample, the 
largest difference in mean colour is A{g — z)q ~ 0.13 (see Ta- 
ble[5]) and for the entire ACSVCS sample, the m ean GC sys- 
tem colour has a range of 0.76 s: (ff-z)o < 1-25 (|Peng et all 
l2006al ). Therefore, stellar populations in GC systems are 
mainly governed by their host galaxy mass (luminosity) and 
environmental effects are less important in determining the 
star formation history of early-type galaxies and their GC 
systems. 

There have been a few previous attempts to detect dif- 
ferences in GC properties depend ing on the environments to 
which their host galaxies belong. ICebhardt fc Kissler-Pati3 
(|l999t ) used 50 mostly early-type galaxies from the 
HST/WFPC2 archive and found no correlation between the 
colour peaks and the host galaxy properties in 15 field sam- 
ple galaxies. However, these researchers found a strong cor- 
relation of GC system peak colours versus galaxy veloc- 
ity dispersion and Mg2 index in their cluster galaxy sam- 
ple. In our field galaxy results, a strong relationship be- 
tween the mean colour and the galaxy luminosity was cer- 
tainly observed (see Figure [9ll. One possible reason why 
ICebhardt fc Kissler-Patia ()l999l ) did not find such a rela- 
tionship in their field galaxy sample is that their galax- 



ies were too luminous (My 



-20) and only small differ- 



ences are to be ex pected for luminous (Mb S, early- 
type galaxies fe.g.. lBaugh et al.ll 19961 : ICole et al.ll2"o0ol ). An- 
other possible reason is that the photometric accuracy of 
the HST/WFPC2 data was not sufficient to detect the 
co rrelation. In fact, the pr e sent work and the research 
of ICebhardt fc Kissler- Patid j 19991 ) have two field galaxies 
in common: NGC 1426 and NGC 3377. When compared 
to iGebhardt fc Kissler- Patid (Il999l). we detected twice as 
many GCs in these two galaxies. Larsen et al.l l|200 also 
briefly investigated the dependency of colour bimodality on 
the environment in an independent sample of 17 early-type 
galaxies obtained from HST/WFPC2, but did not find any 
correlation. As can be seen from our results, the colour bi- 
modality of GC systems mainly depends on the galaxy lu- 
minosity (mass), although this trend disappears at slightly 
higher galaxy lu minosities in low-density regions. 

■Peng et al found a large spread in the specific 

frequency (Sn) of dwarf galaxies from the ACSVCS. Almost 
all dwarfs with a high Sn were locate d within a project ed 
radius of 1 Mpc from M87 (Figure 4 in lPeng et ail (|2008t )). 
whereas no Sn difference between the central region and 



the outskirts of the Virgo cluster was found for intermediate 
luminosity gala xies (i.e. sim i lar to those in our low-density 
sample). While IPeng et al.l (|2008') found a dependency of 
Sn for the dwarfs on the environment, the cumulative GC 
colour distributions of these two groups show no obvious 
difference. It is possible that the effect of environment on the 
colour distribution is so subtle that it could not be detected 
within a cluster environment such as Virgo. 

Regarding the general dependence of galaxy formation 
on environment, many observational and theoretical stud- 
ies have been conducted without including constraints from 
GC system data. It is now well-established that late-type 
galaxies are biased to low-density environments and giant 
ellipt ical galaxies a re preferentially located in galaxy clus- 
ters (|Dresslejll980l ). Cluster galaxies also exhibit lower star 
formation rates than field galaxies at a given r edshift, lu- 
minosity, and bulge-to-disk mass rati o (jBalogh et al. 199i; 
I Lewis et al.l I2OO2I : Gomez et al.l l2003l ) . Early- type galaxies 
in low-density environments also appear to be younger 
(Aage ~ 2 Gyr) and more metal-rich {A[Fe/H] ~ 0.1 
dex) when compa red to their counterparts in dense envi- 



ronmen ts (iKuntschner et al 



2OO2I : iThomas et aLllioOsI ). al- 



though IClemens et al.l ( 2006 ) found no environmental in- 



fluence on the metallicity (but still found younger ages for 
fleld galaxies). These observational findings are contradic- 
tory to the semi-analytic mod el of hierarchical galaxy for- 
mation (e.g. lBaugh et al.ll 19961 . ICole et al, 2000 ), which pre- 
dicts a larger spread of age and metallicity and younger 
ages and lower metallicities (on average) in low-density en- 
vironments, especially in low-lumino sity early-type galaxies 
iBaugh et al.lll996l . lCole et al.ll2000h . 

Our results for the GC systems are s o mewh at in- 
consistent wi t h tho se of IKuntschner et al] (|2002l ) and 
I Thomas et al.l (|2005l ). who found higher metallicities for 
galaxies in field environments ( which is con s istent with 
the semi-analyti c models from iBaugh et al.l 1 19961 ) and 
ICole et al.l (|2000l )). The question remains as to whether GC 
formation history is different from field star formation his- 
tory or does the environ ment have differe nt effects on GCs 
and field star formation. [Harris fc Harrij (|2002) have sug- 
gested that GCs formed somewhat earlier than most field 
stars based on their discovery of a higher Sn at a low 
me tallicity in NGC 5 128. This argument is also supported 
bv iPeng et all (EoO^, who estimated GC formation rates 
from the Millennium Simulation (Sp ringel et al.ii2005i ) and 
showed that cluster formation peaks earlier than field star 
formation. Such results would explain why GCs are more 
metal-poor than field stars in early-type dwarf galaxies. The 
higher Sjv of dwarf ellipticals in the Virgo cluster center 
was successfully explained by the higher peak star formation 
rate (SFR) and star formation surface density (Esfa) in the 
cluster central region when compared to the outskirts and a 
model where GC formation rate CFR oc SFRCSsfr)°'^ ■ 

As described in >i4.4l the colour of GC systems is mainly 
determined by the mass of the host galaxy in both clus- 
ter and field environments. Subtle environmental differences 
also exist; the mean metallicities of GC systems in field 
galaxies are slightly lower than those in cluster galaxies at a 
given host galaxy luminosity (mass), while the fraction of red 
clusters is higher in dense environments. From our findings, 
it is suggested that GC systems in field galaxies form later 
than those in cluster galaxies and/or have a lower metal- 
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licity. Furthermore, we expect that in dense environments, 
cluster formation history is complicated by disturbing or in- 
teracting neighbor galaxies. This is reflected by the larger 
degree of variation in the shapes of the colour d istributions 
that are found in dense environments. While IPeng et al.l 
l|2008l ) found no difference in the Sn of low-luminosity el- 
lipticals (intermediate luminosity according to their termi- 
nology) in their Virgo cluster sample, in the future it would 
be worth comparing the Sn of our sample with that of the 
ACSVCS. 



5.2 Metal-poor GC colour— galaxy luminosity 
relation 

The existence (or non-existence) of a correlation between 
the peak of the blue GC colours and the host galaxy lumi- 
nosity can have implicati ons for different galax y formation 
scenarios. M ajor merger (lAshman fc Zepj|l992l ) and accre- 
tion models l|C6te et aLlligOSl ) would have difhcultv explain- 
ing this correlation because, in the merger model, metal- 
poor GCs in two lower equal mass spiral galaxies are still a 
primary resource of metal-poor GCs in the final ellipticals. 
In the accretion model, metal-poor GCs come from dwarf 
galaxies, so the mean metallicities of GCs in dwarf galaxies 
and metal-poor GCs in ellipticals are more or less the same. 
Many early studies found a strong correlation between the 
colours of metal-rich GCs and the host galaxy luminosity, 
but failed to detect any evidence of a simi lar correlation with 
the mean colours of mctal-poor GCs fe.g..lForbes et al.lll997l : 
Kundu fc W hitrnorc 2001a; F orbes fc Fortell200ll 'l. However. 



Larsen et al.l lj200ll 'l found a weak correlation between the 



mean colours of metal-poor GCs and the galaxy absolute 
magni tude in their HST/WFPC2 sample. iBurgarella et al.1 
ll200ll) also suggested that such a relationship may exist. 
iLotz et al] ([2004l ) found that the slope of the GC peak of 
their dwarf elliptical sample (Afs 1^ —18), when plotted 
against galaxy luminosity, is consistent with that of the 
GC b lue peak in early-type galaxies found bv lLarsen et al.l 
l|200lh . 

Recent studies of the correlation between metal-poor 
G Cs and host galaxy lumi nosity have been c arried out 
bv lStrader et al.l (|2006l l and IPeng et"al] (|2006al l. The for- 
mer compiled blue GC peak data for earl y-type galax- 
ies a n d local spirals from various sources (iLarsen et al.l 



200ll: iKundu fc Whitmor3 l2001al. 12001 bl: iHarrij 1199^ 



Barmbv et aD l2000l : lOlsen et all l2004l l and found a signif 



leant (> 5a) corr elation be tween the blue colour peak and 
galaxy luminositv. IPeng et a l. (2006a) also detected this blue 
peak trend in homogenous ACSVCS data, but it was not as 
strong as that for the red peak. Using their own r elation - 
ship between g ~ z colour and [Fe/H], IPeng et al.l (j2006al l 
found that the slope of the blue peak met allicity- galaxy lu- 
minosity relationship is even steeper than that found in the 
colour-galaxy luminosity plane due to the steeper colour- 
met allicity r«l^^tionsliipin_the metal-poor region (see Fig- 
ure 12 in IPeng et al.l (|2006al )). However, they noted that 
this blue slope varies depending on the adopted colour- 
metallicity relationship (simple stellar population model or 
empirical transfor mation) by a factor of ~ 3. 

In Figure 2 of lBrodie fc Stradeij ()2006l l. the peak metal- 
licities of GC sub-populations are plotted agai nst galaxy lu- 
minosity with the peak positions taken from IStrader et al.l 



(|2004 ) and IStrader et al] (120061 ). In this plot, the peak 
metallicities of both metal-poor and metal-rich GCs appear 
to be correlated with host galaxy luminosity with similar 
slopes. The metal-p oor GC trend can be accounted for with 
an in-situ scenario l|Forbes et al.lll99'i1) in which metal-poor 
GCs formed in the Universe first, and after a sudden trunca- 
tion (possibly by reionization) metal-rich GCs then formed 
along with the bulk of the field stars in galaxies. 

From our results on the blue peak colours alone, it is 
not clear whether or not the blue peak colours are corre- 
lated with the host galaxy luminosity because of the narrow 
galaxy luminosity range in our sample (only ~ 2 mag range 
in Mb compared to a range of ~ 6 mag in the ACSVCS) . To 
observe any trend in the blue peak colour with galaxy lumi- 
nosity, we would need to study GC systems in fainter galax- 
ies {Mb ^ —18, i.e., dwarf elliptical s) in low-d ensity environ- 
ments (similar to those studied bv lLotz et a l. (2004)). This 
would serve to strengthen the correlation between the GC 
blue peak and galaxy luminosity in the higher density Virgo 
and Fornax clusters. It is thus unclear whether the slope 
of the blue (metal-poor) GC peak in field environments is 
different from that in clusters. In fact, Strader ct al. (200^) 
predict that the slope of metal-poor GC metallicity-galaxy 
luminosity plots in dense environments should be steeper 
than that in low-density environments if cosmic reioniza- 
tion caused a quenching of blue GC formation. The role of 
cosmic reionization in GC formation is, however, uncertain. 
Whether cosmic reionization can really truncate GC forma- 
tion and how efficient GC forma tion is prior to the reion- 
ization is still con troversial (e.g. iMuratov fc GnedinI l2010l : 
iGriffen et a"l]|2010l ). 



6 CONCLUSIONS 

• High spatial resolution images of 10 early-type galaxies 
in low-density (field) environments have been obtained us- 
ing the HST/ACS in the F457W and F850LP bands. The 
properties of the GC systems associated with these galaxies 
and their connection with host galaxy properties have been 
investigated. Our results have been compared with those of 
the ACS Virgo Cluster Survey (ACSVCS) in order to study 
the role of the environment in galaxy formation. 

• The GC system properties of our low-density sample 
exhibit trends with respect to the host galaxy luminosity 
that are similar to those in clustered environments. There 
are more total GCs and more red GCs in luminous galaxies, 
while the mean colour and the colour of the red peak GCs 
are strongly correlated with host galaxy luminosity. Colour 
bimodality becomes less clear with a decrease in galaxy lu- 
minosity. 

• The mean colours in low-density (field) regions appear 
to be slightly bluer than those around galaxies of equivalent 
luminosities in high-density regions. The fraction of red GCs 
in field galaxies is found to be lower than that in clustered 
galaxies of similar luminosities. When compared to the trend 
observed in the ACSVCS, the slope of the red peaks against 
host galaxy luminosity is steeper. In other words, colour 
bimodality disappears more quickly as galaxies get fainter, 
in low density environments. 

• Luminosity functions of the GCs in most of our sample 
galaxies are well-fit by a Gaussian function whose fitting 
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parameters (turn-over magnitude and dispersion) agree well 
with those of the ACSVCS. Wc find no independent evidence 
for a trend of dispersion in the fitted Gaussian luminosity 
functions against galaxy luminosity. 

• We have investigated the possible origins of the unex- 
pected excess of GCs in NGC 3377 near g 25.5, which 
deviate from a normal Gaussian luminosity function. The 
origin of this population is still not fully understood, but 
some of clusters are likely to be background galaxies from 
visual inspection and some arc possibly diffuse star clusters. 

• At a given luminosity, the mean colours of our GC sys- 
tems in low-density environments are slightly bluer {A{g — 
z)o ~ 0.05) than their counterparts in high-density envi- 
ronments. By assuming various combinations of age and 
employing simple stellar population models, this colour off- 
set corresponds to a metallicity difference of A[Fe/H] ~ 
0.10 — 0.15 or an age difference of at least Aage ~ 2 Gyr on 
average. Therefore, GCs in field galaxies appear to be either 
less metal-rich or younger than those in cluster galaxies of 
the same luminosity. 

• Whilst a correlation between the blue peak colour of 
GCs and host galaxy luminosity was found in the Virgo 
cluster (ACSVCS), no such correlation was detected in our 
low-density galaxy sample. It is not clear whether the ab- 
sence of such a trend is because our galaxy luminosity range 
is simply too narrow or because intrinsically the slope of the 
relationship depends on the environment. For a more defini- 
tive determination, observations of GCs in dwarf ellipticals 
in low-density environments are needed. 

• The greater variation in the shapes of the colour distri- 
butions for GC systems in the Virgo cluster sample could 
imply that more complex galajcy formation processes (e.g., 
interactions/harassment with adjacent galaxies) arc taking 
place in the galaxy clusters. The higher fraction of red GCs 
in cluster galaxies also supports this possibility. 

• Although we found that the galaxy environment has 
a subtle effect on the formation and metal enrichment of 
GC systems, host galaxy mass is the primary factor that 
determines the stellar populations of both the GCs and the 
galaxy itself. The processes which determine GC formation 
must therefore bo largely common across a range of galaxy 
environments. This remains a challenge for GC formation 
theories. 
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